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Abstract 

We compare two opposite ways of performing a 3D reduction of the 
two-fermion Bethe-Salpeter equation beyond the instantaneous approxi- 
mation and Salpeter's equation. The more usual method consists in per- 
forming an expansion around an instantaneous approximation of the prod- 
uct of the free propagators (propagator-approximated reduction). The 
second method starts with an instantaneous approximation of the Bethe- 
Salpeter kernel (kernel-approximated reduction). In both reductions the 
final 3D potential can be obtained by following simple modified Feynman 
rules. The kernel-approximated reduction, however, does not give the cor- 
rect scattering amplitudes, and must thus be limited to the computation 
of bound states. 

Our 3D reduction of the three-fermion Bethe-Salpeter equation is in- 
spired by these results. We expand this equation around positive-energy 
instantaneous approximations of the three two-body kernels, but these 
starting approximations are given by propagator-approximated reductions 
at the two-body level. The three-fermion Bethe-Salpeter equation is first 
transformed into a set of three coupled equations for three wave functions 
depending each on one two-fermion total energy, then into a set of three 
3D equations and finally into a single 3D equation. This last equation 
is rather complicated, as it combines three series expansions, but we use 
it to write a manageable expression of the first-order corrections to the 
energy spectrum. 
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PACS ll.lO.Qr Relativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 
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1 Introduction 

The Bethe-Salpeter equation |jl|, || is the usual tool for the study of relativis- 
tic bound states. The principal difficulty in the treatment of this equation 
comes from the existence of unphysical relative time variables. In the two-body 
Bethe-Salpeter equation the simplest way of eliminating the (in this case single) 
relative time variable consists in performing an instantaneous approximation 
of the Bethe-Salpeter kernel (kernel-approximated reduction) : one replaces the 
complicated kernel given by field theory by a simpler expression, independent 
of the initial and final relative energies in the momentum space representation. 
A 3D equation (Salpeter's equation |^ ) follows immediately. It is also possible 
to get Salpeter's equation from an instantaneous approximation of the product 
of the two free propagators: one replaces this product by a 3D propagator times 
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a delta (the " constraint" ) fixing the relative energy (propagator-approximated 
reduction). 

The derivation of Salpeter's equation using a kernel-approximated reduction 
appears to be the most natural one. For the calculation of the higher-order 
correction terms (towards a theoretical equivalence with the original Bethe- 
Salpeter equation, at least for what concerns the bound state energy spectrum), 
the propagator-approximated reduction has been preferred by most authors [4- 
19] including us 18, |l^. We found an example of a kernel-approximated 



reduction in articles by Phillips and Wallace |gO| , who tune their instantaneous 
approximation of the kernel in order to make the corrections vanish up to a 
given order. It is however our tentative of 3D reducing the three- fermion Bethe- 
Salpeter equation pl| which led us to consider this alternative approach (see 
below). 

This work is divided into two parts. In the first part we present the two 3D 
reduction methods of the two-fermion Bethe-Salpeter equation. We do not insist 
on the propagator-approximation reduction as we already closely examinated 
it elsewhere (l7[ |l|, |l|l, and we focus on the kernel-approximated reduction. 



comparing it with the other method. 

We meet a first difficulty coming from the fact that the kernel-approximated 
reduction does not lead directly to an symmetric potential. In the case of a 
total-energy depending potential like here, this does not necessarily imply that 
the final energy spectrum will not be real. It is however interesting to get a 3D 
equation with a symmetric potential by performing a supplementary expansion. 
The resulting double expansion (triple if we consider also the higher-order irre- 
ducible graphs contributing to the kernel) does not complicate the enumeration 
of the contributions. It is indeed possible, as in the propagator-approximated 
reduction, to get the final potential by slightly modifying the Feynman rules 
used in the calculation of the 4D off mass shell transition operator. 

The most concise way of presenting the two reduction methods consists in 
writing their respective 3D off mass shell transition operators in terms of the cor- 
responding 4D off mass shell transition operator of field theory. Two reductions 
which share the same 3D off mass shell transition operator correspond indeed 
to 3D equations which can be transformed into each other at the 3D level. In 
the propagator-approximated reduction, the 3D operator is obtained from the 
4D one by applying a constraint which is also satisfied by the free solutions, so 
that the corresponding 3D equation gives not only the bound states, but also 
the correct physical scattering amplitudes, although the original homogeneous 
Bethe-Salpeter equation was only valid for bound states. We loose this bonus 
with the kernel-approximated reduction. 

We close this first part with the calculation of the second-order contributions 
to the 3D potential using the kernel-approximated reduction method, in the case 
of a simple kernel given by a one-photon exchange graph (in Feynman gauge). 
Each term is given by the residues of two poles in the initial and final relative- 
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energy complex planes. The poles can come from the propagators only, or one 
of them can come from the kernel. The contributions of the kernel poles are not 
negligeable, at least in this case. Beyond the ladder approximation, however, 
they could be partially cancelled by the contributions of the first crossed term. 

In a recent work [ pT| , we proposed a three-fermion Bethe-Salpeter equation with 
three instantaneous, cluster-energy independent and positive-energy kernels. 
Such an equation can easily be transformed, without supplementary approx- 
imation or perturbation expansion, into a 3D equation. The instantaneous two- 
body kernels were chosen in order to reproduce the two-body potentials obtained 
with a propagator-approximated reduction in the pure two-body case. With this 
choice, our 3D equation satisfies a cluster separability requirement: "switch- 
ing off' two of the three mutual interactions leads to a free equation for one 
fermion, plus a 3D equation for the two other fermions, equivalent to the exact 
Bethe-Salpeter equation for these two fermions. The remaining discrepancy (at 
least in the predicted energy spectrum) between the exact three-fermion Bethe- 
Salpeter equation and our three-fermion 3D equation must then be considered 
as a three-body relativistic effect, the computation of which is the purpose of 
the second part of the present work. Around a positive-energy instantaneous 
approximation of the three two-body kernels (the three-body irreducible terms 
being entierely included in the perturbation) , we perform three consecutive ex- 
pansions of the correction terms: the three-fermion Bethe-Salpeter equation is 
first transformed into a set of three coupled equations for three wave functions 
depending each on one two-fermion total energy, then into a set of three 3D 
equations and finally into a single 3D equation. The resulting expression for 
the 3D equation is rather complicated, as it combines three series expansions, 
but we use it to write a manageable expression of the first-order corrections to 
the energy spectrum. The contributions of the poles of the propagators to this 
energy shift are easy to compute in terms of the Bethe-Salpeter kernels with 
two of the three fermions on their mass shell, in the spirit of Gross' spectator 
model |2^, |23| . The contributions of the singularities of the kernels themselves 
(which could be important) remain to be computed but this should be done in 
the framework of more specific models. 

What about cluster separability? We should not be surprised by uncor- 
rect cluster-separated limits, as the three-fermion homogeneous Bethe-Salpeter 
equation is valid only for three-body bound states. Correct limits would again 
be a bonus. The cluster-separated limits of our 3D equation are two-fermion 3D 
equations built with the potentials corresponding to the starting instantaneous 
two-body kernels. They are thus correct if these potentials are the full poten- 
tials corresponding to reductions of the two-fermion Bethe-Salpeter equation. 
Another choice, like that based on the simpler Born approximation, would lead 
to a 3D equation whith uncorrcct cluster separated limits, but could still be 
used in the computation of three-fermion bound states. 



4 



2 The two-fermion problem. 
2.1 Notations. 

We shall write the Bethe-Salpeter equation for the bound states of two fermions 
i as 

$ = G"X$, (1) 

where $ is the Bethe-Salpeter amplitude, function of the positions xi,X2 or of 
the momenta pi,p2 of the fermions (according to the chosen representation). 
The operator K is the Bethe-Salpeter kernel, given in a non-local momentum 
representation by the sum of the irreducible two-fermion Feynman graphs. The 
operator G° is the free propagator, given by the product G^G" of the two 
individual fermion propagators 



G° = r^— . A (2) 



where the hi are the Dirac free hamiltonians of eigenstates itE'i : 

hi ^ di .pi + f3imi (i = l,2). (3) 

E, = .fh^=if,+m^)i. (4) 



We shall define the total (or external, CM, global) and relative (or internal) 
variables: 

X^-{xi+X2), P=pi+P2, (5) 

X = Xi-X2, P^^{Pl-P2)- (6) 

and give a name to the corresponding combinations of the free hamiltonians: 

S = hi+h2, S=l(/li-/l2). (7) 

2.2 3D reduction using an instantaneous approximation 
of the propagator. 

The free propagator G° will be approached by a carefully chosen expression G*, 
combining a constraint fixing the relative energy, and a global 3D propagator. 
There exists an infinity of possible combinations [4-19] which would be equiv- 
alent if all correction terms could be computed. Our choice in this work will 
be 

G'=A++ J dpoG^ipo) - -2t7rSipo-s)A++g°P,l32. (8) 
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where A"*"+ is the free positive-energy projector: 

A'^^AlAi Af = ^%^- (9) 
and 50 the 3D propagator: 

This choice of is not necessarily the best one for working in the pure two- 
fermion framework, but it is the simplest one for the presentation of our two 
and three- fermion results. We shall write the free propagator as the sum of the 
approached propagator, plus a remainder: 

G" = G'" + G". (11) 

The Bethe-Salpeter equation becomes then the inhomogeneous equation 

<^> ^G^K^ = {G^ + G^)K<i> = + G'^K<S>, (12) 

with 

* = G^K^. (13) 

Solving (formally) the inhomogeneous equation (|lj) and putting the result into 
(H, we get 

* = G'^/^(l - G^K)-^^/ = G*A'^* (14) 

where 

= K{1 - G^K)-^ = K + KG^K + ... (1 - KG^Y^K (15) 



= K + KG^K^ = K + K^G^K. (16) 



obeys 



The reduction series (|15|) re-introduces in fact the reducible graphs into the 
Bethe-Salpeter kernel, but with G° replaced by G^. Equation (|lj) has the 
same bound state spectrum as the original Bethe-Salpeter equation. 

In the inversion of [X—GrK) and similar operators below, we have to assume 
that eq. {l — G[iK)^ — Q has no solution for the values of Pq considered (i.e. 
the values belonging to the bound state energy spectrum of the Bethe-Salpeter 
equation). 

The relative energy dependence of eq. (|lj) can be easily eliminated: 

* = %o-s)V' (17) 

and ip obeys: 

V' = 3o^^V' (18) 
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where is proportional to the positive-energy part of with the initial 
and final relative energies fixed to s : 

V'^ = -2z7r A++ f3if32K^{s, s) A++. (19) 

In less compact but more precise notations: 

pip2K''{s,s) = J dp'odpoS{p'o-s)l3il32K^{p'o,po)S{po-s). (20) 

Note that we write {p'q,Po) but (s,s), as we keep s in operator form. 

These 3D reductions can also be described in terms of transition operators. 
The 4D transition operator is 

T = K + KG°K + --- (21) 

and can be obtained by keeping only the part of G° in it. We have 
also 

T = Kil-G^Ky^ = K{l-G^K-G^K)-^ = K{1-G^K)-'^{1-G^K{1-G^K)-'^)-'^ 

= K^{1-G^K'^)-^ = + K^G^K^ + ■■■ (22) 
so that the 3D transition operator is given by 

+ g°V^ + ■■■= -2i7r/3i/32 (23) 

with 

= K^++{s, s) + K^++{s, s) K^++{s, s) + ■ ■ ■ = T++{s, s) (24) 

K^++is,s) = /3i/?2A++/3i/32i^^(s,s)A++ = ^ p^p^V^ (25) 

T++{s, s) = Pi(i2 A++ Pi(i2 T{s, s) A++. (26) 

We see that the 3D transition operator is a constrained form of that of field 
theory. Both operators become equal to the physical scattering amplitude when 
both fermions are on their positive-energy mass shells. Conversely: 

T = K'^ + K^G^K^ + K'^G^K'^ + K'^G^K^ + ■ ■ ■)G^K'^ 

= K^ + K^G^K^ + K^G^T^G^K^. (27) 

A Bethe-Salpeter equation leading directly to the same 3D reduction can be 
obtained by replacing the kernel K by the instantaneous positive-energy ker- 
nel K'^^~^{s, s). In this case the 4D transition operator T becomes equal to 

T++{s,s)=T(^. 
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The first-order correction term to the energy is 

Po - ^ -2iTr < A++ /3i/32 {KG^K){s, s) A++ > . (28) 

We must not forget that P^^\ {KG"'K){s,s) and the zero-order wave func- 
tions are in general total energy dependent, so that ( |2^ ) is in fact a numerical 
equation in Pq (or a matricial equation in case of degenerescence) . If the Pq 
dependence of the zero-order potential is of higher-order, one has only to take 
the Po dependence of Pq°^ into account. 

2.3 3D reduction using an instantaneous approximation 
of the kernel. 

We shall now consider an approximation of the Bcthe-Salpeter kernel: 

K ^ K° + K^. (29) 
The Bethe-Salpeter equation becomes 

$ ^ G^K^^ + G°K^<^ (30) 
$ = (1 - G^K^T^G^K^'^ = G^K°<!> (31) 

with 

- G" + G°i^«G" + G"K^G"K''G'> + • • • ee G" + G^^ . (32) 

If we now specialize to an instantaneous positive-energy kernel {K^ inde- 
pendent of Po and equal to /3i/32A++/3i/32i^°A++ ), eq. (^) leads to the 3D 
equation 

= (s"+5'^^)l^0 (33) 

with 

4> = K++ j dpo $(po), V = -2z^/3i/32i^°, (34) 

= A++ / dp'o dpo G^«(p^,po) /3i/32 A++. (35) 
The first-order energy shift is, with the replacement g^^V(t) ~ g^^ {Pq — S)(f) : 

Po-p!,''' = < iP"-S) J dp',dpoG"{p'o)K^{p'o,Po)G"{po)Pip2{P'~S) > . 

(36) 

Here again, all elements of ( |36| ) are total energy-dependent. Although ( p6| ) 
is symmetric, the perturbation potential {Pq — S)g^^V of (^) is not. This 
does not exclude the possibility of a real energy spectrum: if we admit an 
energy dependence of the potential (which we can not avoid), we have also 
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the possibility of performing an infinity of rearrangements of the equation. We 
can for example render the potential symmetric by expanding eq. ( p3| ) in the 
same way as eq. ( p^ of section (2.2). We can equivalently start with the 
corresponding 4D equation (pi]): 



$ = (G° + G^^)X°$ = x + G^^-ftr°«' (37) 

X = G"/v"<i> = G°i^"(l - G^^ i^°)-i X = G^K^'x (38) 

where 

=K'^ + K°G^'^K'^ + • • • (39) 

is a symmetric positive-energy instantaneous potential, since it begins and ends 
with K'^ . The 3D equation is then 

77 = g°y^r;, 1/^ = -2z7r/3i/32if^, V=[dpoxiPo) (40) 



and the first-order energy shift is still given by (^6|). As for the propagator- 
approximated reduction, we can represent the various contributions by Feynman 
graphs. We can indeed get by considering the transition operator 

T = {K° + K'^) + (K" + K^)G°{K° + K^) + ■■■ (41) 

and skipping some terms according to the two rules: 

- Each term must begin and end with K'^ 

- There must be at least one between two consecutive K'^. 
The 3D transition operator is now —2i7T(3i/32T^ with 

= A'^ + X^G°/^^ + • • • = i^^(l - G°if^)-i 

= K"il - G^^i^")-\l - G°i^°(l - G'^'K^rY' 
= K°il - G'^'^K^ - G°K°)-^ = K"il - G'^kY^ 
= i^"(l-(l-G"i^^)-iG"i^")-i 
= K%1 ~ G^K^ - G°K")-\l - G"k" ) = K"{1 - G°K)-\l - G^K'^ ) 

K°{l + G^T){l-G^K'^) (42) 
Replacing then by K~K° and r(l-G''A:) by K gives 

T'^ ^ K" + K"G"K" + K"G" T G"K°. (43) 

In contrast with what happens in the propagator-approximated reductions, the 
on mass shell restriction of this operator is not the physical scattering amplitude 
of field theory. In fact, the Bethe-Salpetcr equation (|l]), without inhomogeneous 
term, is a bound state equation, and for the bound states we are only interested 
in the poles of (H3). These poles are the poles of T, unless the integrations on 
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p'q,Po make some residues vanish. Getting also the correct physical scattering 
amplitudes was a welcome bonus of the propagator-approximated reductions. 
If we want T in terms of , we can write ( p2| ) as 

T^(l -G°i^^ )-i = i^"(l + G°T) {K - K^ ){1 + G"T) 

= T - K" {1 + G^T) = [l- K^G° )T - (44) 
T = (1 - K'^ G" T^^' (1 - G°i^^ + (1 - if ^ G" ii^^ 

= (1 + r^G°) (1 + G°T^) + (45) 

with 

^ kR + K^G^K^ + ■■■. (46) 
We have a large freedom in the choice of K'^. We could choose 

K'> = /3i/32 A++ /3i/32 if(s, s) A++. (47) 

We could also keep only the ladder term of K, or a part of it in (^7|), like the 
Coulomb part of the one-photon exchange contribution. Wc could use (in 
this case we expect that the remainder of docs not finally contribute to the 
total energy spectrum) or an approximation of it (one and two photon exchange 
terms, for example). Phillips and Wallace |Q suggested to choose K'^ in order 
to annihilate exactly or up to a given order. If we choose to annihilate it 

at first-order, the first-order energy shift ( |3^ ) vanishes and we get 

K" = (^/3i/32 {gT'^^^ J dp'„dpoG''{p'o)K{p'„,po)G^p,) (3,(32 {gT'^^^- 

(48) 

The skipping of terms in ( ^H ) is then improved, as we must now have at last 
two between two consecutive if°. If we choose to annihilate 5^^^ exactly, 
we get, using (Isl), ^ and the relation G^^ = G°T^G° : 



^/3i/32 {gT'^^^J dp',dpoG"{p'„)Tip'„,po)G'ipo)^)i^)2 (5°)-'A++ = 

(49) 

which is the A++ projection of the expression given in (they use Salpeter's 
propagator and write thus A++— A^^ instead of A++). In this case, gives 
the correct physical scattering amplitudes when P'^ ^ E1+E2, as the operators 
((7°)~^ kill the contributions of the singularities other than the poles of G". 

2.4 One-photon exchange. 

Let us now compute ( p8| ) with K given by the one-photon exchange graph, in 
Feynman gauge: 

^ 2^e^ 71-72 , s 
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Writing p'—p = k and leaving out trivial factors, we have to compute the integral 



/ 



dpio dko 



1 



Pw 



■ le 



Po - Pw -ko - E'2+ie 
1 



k^ + ie Pw - Ei+ie Pq - pio - £'2 + ie 



(51) 



Performing the integration in pio by closing the integration path clockwise, we 
get the residues of the poles at pio = E[ — ko, Ei : 



I = -2iTT 



dko = 

(ri) k^- k^ + ie 



1 



Po-E[-E2+ko+' 



+ 



Po-E[-E'2+ ie E[-Ei-ko 



1 



El- E[+ko + ie Pq - Ei - E'2 - ko + ie 



1 



(52) 

Pa-Ei-E2+ie j ^ ' 

The possibility of kg being equal to E[—Ei, leading to a double pole, has been 
taken into account by leaving an interval of length rj around this value in order 
to take the 77 — » limit afterwards. This can also be done by replacing the 
poles ±(E[^Ei — ko±ie)^^ by their principal parts, or equivalently by adding 
a iii5{ko — E[ + El) to each one. For the integral in fco we shall close the 
integration path counterclockwise in the first term and clockwise in the second 
one. In each case we shall have a pole of the propagator at ko = E[—Ei (the 6 
reduces the contribution of this pole by a factor 2) and a pole of the kernel at 
A;o = — 1^1 and fco = |fc| respectively. The final result is: 



I = -47r^ 



1 



1 



Po-E[- E'2 (^E'l - Elf -P Po-Ei- E2 




El - \k\ )(Po -El- E2) {E[ -Ei + \k\ )(Po - E[ 



1 + 



E'2) 



Po-E[-E2- \k\ Po-Ei-E'2- \k\ 

(53) 

Here again, we must take the principal parts of the poles in \k\ = \E[—Ei \ in the 
future integrals in k. The first term comes from the residues of the poles of the 
propagators only, while the two other terms combine the residues of one pole 
of the propagators and one pole of the kernel. Closing the integration paths 
in the opposite directions would interchange the two fermions in the first term 
and modify the two other terms accordingly. At lowest order, the expression 
between brackets becomes 



1 



1 + (2Po 



E'l — E'2 — El — E2). 



(54) 
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We could define K'^ as the kernel which leads to (p3|), to the first term of ([5^), or 
to the first term of ( |5^ ) with {E'l—Eif' —}? replaced by simply —k^ (Coulomb 
potential). We see anyway that the relative time dependence of the kernel 
contibutes more to (|5^) via the residues of its poles than directly. We should 
however consider also the contribution of the lowest-order crossed graph, which 
could partially cancel (^). We know indeed that, in well chosen propagator- 
approximated reductions, all higher-order contributions cancel mutually at the 
one-body hmit g |l|. 

3 The three-fermion problem. 

3.1 Three-fermion Bethe-Salpeter equation. 

The Bethe-Salpeter equation is now: 

$ = [g;G°A'i2 + G^Ggif 23 + GlGlK^i + G;G§G°A'i23] $ (55) 

where Ki2'i is the sum of the purely three-body irreducible contributions. Let 
us first consider the case of three two-body positive-energy, instantaneous and 
cluster-energy independent interactions K?^^ (i.e. containing A^'*' projectors 
and independent of p^^Q^pij^^Pijo). The Bethe-Salpeter equation becomes in 
this case 

$ = G?G°G° /3i/32/?3 [ Vi2 ^12 + 1^23 ^-23 + 1^31 V'31 ] (56) 

Zm 

V,, = ~2i^ PSjKl , {Pm) = Pk J dp„o (57) 

This leads to a set of three coupled integral equations in the tjjij . We shall search 
for solutions analytical in the Im(p/jo) <0 half planes and close the integration 
paths clockwise in these planes. The only singularities will then be the poles of 
the free propagators. Performing the integrations with respect to the pijo gives 
then 

A++ 

fpuiPso) = 775 7^ , ^ , . [ Vi2'012(P3o) + ^23^(^1) + 1^31 V'31 (^2) ] 

[Pq - S) - (p3o - hs) +ie 

(58) 

and similarly for -023 and '031 (from now on, we shall often write the equations 
for the (12) pair only, understanding the two other equations obtained by cir- 
cular permutations on the ijk indexes). Solving (js^ ) with respect to V'i2(P3o) 
gives a wave function which is analytical in the Im(p3o) < half plane, and 
confirms thus the existence of such solutions. Furthermore, equation ([58| ) shows 
that the three projections A^V-'y (^fe) are equal (let us call them ip ) and satisfy 
the 3D equation 

= - -[V,2 + V23 + V,^]^. (59) 

/t) — £/l — £/2 — -C'3 
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Moreover, ip is the integral of the Bethe-Salpeter amphtude with respect to the 
relative times; 



■0 = TT— dpo^ = — dpiodp2odp30^{PlO+P20+P30- Po)^- (60) 



in order to recover, at the three 2+1 separated clusters limits, the propagator- 
approximated 3D reductions (subsection 2.2 above). 

For the general three-fermion Bethe-Salpeter equation ( p5| ) it seemed to be 
a good idea to try an expansion around a set of positive-energy instantaneous 
approximations of the two-body kernels. We first tried to start with approxi- 
mations giving the 3D equation (E^ in a more direct way than above, like 



correction terms spoil the result: with (|6^) or (|6^) they contain supplementary 
zero-order contributions. With ( |64[ ) they contain no supplementary zero-order 
contributions, but the correction terms can not easily be sorted by increasing 
order. We must thus proceed more carefully. We finally decided to start with a 
positive-energy, instantaneous and cluster-energy independent approximations 
of the two-body kernels, like (^) or its Born approximation, adding the third 
projector in order to avoid unnecessary complications, and to proceed in 
three steps, inspired by the demonstration ( p5| ) -(p9| ): 

i) Transform the Bethe-Salpeter equation (BSFinto a set of three coupled 
equations for the V'ij(Pfco)- 

ii) Transform this set into a set of three coupled 3D equations for the tjjij{hk) 
(these are not the Faddeev equations). 

iii) Transform this last set into a single 3D equation for V'j defined as the 
mean of the ipiji^k)- 

3.2 Three coupled equations for three cluster-energy de- 
pending amplitudes. 

Let us first write the Bethe-Salpeter equation ( ^5|) in a more compact form: 





In we proposed the choice 



with 




(65) 



with 




(66) 
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K = {Gl)-'Ki2 + {G'i)-'K23 + iG'^r'K3i + Kus- (67) 

With these notations, we can formaUy perform the same manipulations as in 
the two-fermion problem (subsection 2.3): 

K = K° + (68) 

$ = G°i^"<I> + G°A'^$ (69) 

<^^{1-G°K"Y^G^K°^ = G'^ K°<^ (70) 

with 

= G" + G°if«G" + G'^K^G'^K'^'G'' + • • • = G" + G^^. (71) 

bmce eq. @, ends with ^0$, we can transform it into a set of three coupled 
equations for the ipij{pko)j using the definition (p7|): 

■012(^30) " X! / ^PkoGi2,zj{P30'Pko) K°j i:,j{pko) (72) 



G^2,Mo,Pko) = P3 [G°Mo)] ' I dp'^^odp^^oG'^ip'o^po) Pk (73) 



with 



3.3 Three coupled 3D equations. 

Let us now compute the contribution of G" to ([73|). We get 

— 2i7r 

G?2.12(P30,P30) = ^^^1 /^l/^2 ^(P30 -P30) (74) 

G;2,23(P30,Pio) = /33G;(pio)G!](Po -Pio -P3o)/3i (75) 

G?2,3l(P30,P20) = /33G°(p2o)G;(Po -P20 -P30)/?2 (76) 

The product of propagators ( |75| ) will be included in a clockwise path integral 
in pio, to which only the pole of G]'(pio) at pio = hi contributes. This is no 
more true in the general case, but we can isolate a part of G5(pio) which leads 
to the same result: 

GUpio) = Gf(pio) + Gf (pio) (77) 



with 



so that 



Giipw) = -217: 5ipw - hi) Pi (78) 

Gf (Pio) = \ Pi (79) 

PiQ - hi - te 

G'uaMo^Pio) = jTS u,^- P2P3S{pio-hi) (80) 
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and similarly 

2i7r 

G'i2,3i(P30'P2o) = 7^ — -r-, j-r—- PzPi S{P20 - ^2)- (81) 

(^0 - j) - (P30 ~ "3j + if 

Putting all together, we have finally 
— 2?7r 

[/32/33 if2°3'/'23(ftl) + /33/3iif3"iV3l(/^2)] 



+ f dpw Gf (pio) G^(Po - Pio - P30) Pi ^23 ^-23(^10) 



+ J dp2o /33 (p2o) G;(Po - P20 - P30) P2 Kl, V31(P20) 

+ E / ^^"^-0 ^^"2!!, bso: Pfco) (Pko) (82) 

If we neglect the last term of (^2|) the two preceding terms vanish also by 
integration, and we recover the result of the positive-energy instantaneous ap- 
proximation. We shall thus move the first term of (|8^) to the left-hand side, 
consider the two next terms (second line) as the principal contributions, and 
the three last terms as perturbations. The difference between this approach 
and that based on ( |6^ lies in the fact that we do not replace immediately p'^Q 
by in the first term. Wc get 

■012(P3o) = .912(^30) [V23'023(/il) + 1^31 ■031 (/i2)] 

+ 512(^30) [ffl2(P3o)]"^ / (^PkoGi2^,j{p'3o,Pko)K^ji',j{pko) (83) 

k 

where G^^^ is defined from G^^^ by including the 3th and the 4th lines of 
( p2[ ) in the last one, and 

512 (Pso) = TF, — — X / , — r~r~~ (^^) 

[Pq- S - V12) - (P30 - /J3) + «e 

5?2(P3o) = T5 — ?^ — A — r~r~~ (^^) 

The iterations of ( p3| ) lead to 

-012(^30) = 512(^30) [^23-023(^1) + ^31-031(^2)] 

+ 9i2ip'so)[9i2iP3o)r' E / dpkoG^,';,^ip',^,pko)K^j 
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gijiPko) [ Vjk tpjkihi) + Vkt ipki{hj) ] . (86) 

defining G^'^ as the sum of tlie iterations of G^^ . 

We can now take (|6|) at P30 = /i3, multiply at left by [5i2(/i3)]^^, and move 
the term in Vi2'0i2(^3) back to the right-hand side: 



Pn-S 



Vi2?Al2(^3) + ^63^(^1) + V3iV'3l(/l2) 



with 



+ ^ Mi2,i-j [Vjkil^jk{ht) + Vktipki{hj)] 

k 

Mi2,rj = / dpknGl^^Ah3,pka) K^j gi]{pkn)- 



(87) 
(88) 



3.4 One 3D equation. 



Let us write ( p7| ) and the two other equations obtained by circular permutations 
in matricial form. We shall define 



^ij{hk), 



i n = 1 - n. 



(89) 



The matrix 11 is the projector (n'^ = n) on the u'^ — \ vector (of norm 3). We 
have 

n* = -(/^u, V = - (*^ + + ^^)- (90) 
3 



We can write eq. (|8^) in the form 

The action of the projector 11 on A'P is 

n 



Pq-S 



V23, etc.. 



hm" = Af", A/r"n = 



p^-s 



V, V = Vi2 + 1^5 



23 



31- 



(91) 



(92) 



We shall now transform the matricial equation (|91|) into a scalar equation for 
if). The expansion of (|9l| ) with respect to 11^' gives 

n * = n Af ( 1 + TiAf + TlAf IlAf + • • •) n ^' (93) 

in which HA/ — liM^. Taking the scalar product with u leads finally to 



— 5— V + - u" [M^ + M[IiM" + IIM^ liM^ + •••]) u 

-Po ^ S 3 



V-. (94) 



In this equation, V contains the projector A"'"++ at left and at right, while 
Af ^ contains it at right. We can thus trivially transform (^) into an equation 

for A+++V-'. 
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3.5 First-order energy shift. 

The 3D reduction above is rather comphcated, since it impHes three levels of 
series expansions. The computation of the first-order energy shift remains how- 
ever tractable. The energy shift is, keeping only the terms which contribute at 
first-order: 

Po - P^"^ =<iPo~S)\u^{l + Af°n ) Af^ u > . (95) 
We can make the replacements 

1 + M°Tl = 1 + V M° (96) 

Po — S 

as the contributions of the first and third terms cancel mutually. We remain 
with 

Po - Pq'"^ = < (Po - S*) i M° M^u > 

k' M 

The contribution of the terms with i' j' = 12 is then 

Ai2 = < ^12 E / '^P^^ (^3' Pfco) 9^, (Pko) {V-V,,)> (98) 

k 

At first order, we can replace Gj^^j by the non-iterated G^^-^j and further by 
Gi2%, as the 3th and 4th lines of (|8^) will lead to null integrals. Using the 
definitions ([7l|), ( [73[ ) we get 

Ai2 = < ^12 ^ / S{p',o - hM dpo /33 [G°Mo)] 

k 

G^''(Po,Po) Pk g^JiPko) iV-V,,)> (99) 



Replacing V by Pq-S and G^^'^ by G°K"G° leads to 

Po-p!,"^ = ^< I dp'„dpoVF{p'o)G''{p'o)K^{p'o,po)G°{po) (3i(32(}3Vi{po) > 

(100) 

with 

ViiPo) = J2 iPko) g-,\hk) (101) 

k 

VAp',) = -2tn J2 V^'f ^P'k'o - hk>)(3k' [G^'(Pfc'o)] (102) 
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In writing (102), we mean that [G°, (Pfc/g)] ^ cancels the G\, (Pt-'n ) con tained in 



G°(Po) before the replacement of p'f.,Q by hk' ■ The energy shift (100) contains 
36 terms: there are 4 terms in (if we do not forget the irreducible three- 
body term) times 3 terms in Vi and 3 terms in Vf ■ Our reason for introducing 
these Vi and Vp is the fact that they both could be replaced by V if 
were instantaneous. 

The singularities to consider when performing the integrations are the poles 
of the propagators, the singularities of and the poles of V/. It is always 
possible to choose the integration paths, in each term, in order to avoid contri- 
butions from the poles of Vi . The contributions of the poles of the propagators 
are then easy to compute, as gij{pka) in ( |101| ) is to be taken at pko = hk- We 
get 

plus the similar R23 and R31 and a contribution R123 of the irreducible three- 
body kernel. By V-^'"^ we denote — 2z7r/3i/32i^f| with the initial fermion a, the 
final fermion b and the spectator fermion 3 on their positive-energy mass shells. 
We have to consider 9 terms in 4 different groups: 



fc' = fc = 3 

k' = lor 2, fc = 3 

fc' = 3, fc = lor2 

fc' = lor2, fc = lor2 



ba ^ 11, 12, 21 or 22 
ba = fc'l or fc'2 
ba = Ik or 2 fc 
ba = k'k 



We still have some freedom in the choice of ba which we could use to get 
a real energy shift, for example by choosing a = 6 in the three first groups. 
The contributions of an instantaneous part of if/l would sum into a 

< V-^^ > but we could also include this directly in Ki2- The corrections 

corresponding to the differences between the various V^^°' are a step towards 
Gross' spectator model |2|, |2^ . 

The principal task remains however the calculation of the contributions of 
the singularities of the Kjj, for the integration paths defined by our choices of 
ab above. As we saw in section 2.4, these contributions could be as important 
as these coming directly from the relative time dependence of these kernels 
(they could however be partially cancelled by the contributions of two-body 
non-ladder terms and three-body terms jlj, ^ ^ ). A general expression 
would be obtained by removing the delta part of the propagators containing 
the poles considered in (103), using ([73|). To go further, however, we should 
specialize to specific three-fermion problems. It must be noted that the ratio 
dijiPko) dij^ihk) of (101) is no more 1 in the terms containing G^{pko)- If 
the singularities of in the complex pko plane lie far from the positive- 
energy mass shell, this ratio can be approximated by —{pkQ — hk)^^g^j^{hk) 
and takes a relatively low value, suggesting us to neglect the contribution of the 
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corresponding G^{pko) to ( |100[ ). We could then make the replacement 

Viipo) ^ = ~2t7r f^k [GUpw)] '^5{Pko - hu) + (104) 

k^u k—u 

in ( |100| ). We isolated the contributions of the unconnected terms k — u (like 
fe = 3 combined with A;' = 3 and ) ■ The resulting energy shift is now 
manifestly real. 



3.6 Cluster separability. 

Physically, if we "switch ofP' the (23), (31) and (123) interactions, we are left 
with a pair of interacting fermions (12) and a free fermion (3). We would like 
to read this on our final 3D equation, however: 

- The cluster-separability condition may be satisfied by an exact 3D reduc- 
tion, but not by its truncated forms. 

- Even with an exact 3D reduction, the cluster separability property is a wel- 
come bonus, not a requirement: the Bethe-Salpeter equation (|55|), from which 
we started, is indeed correct for three-fermion bound states only. 

In our calculations above, the cluster-separability property is in general 
spoiled by our iteration of (p3|), which starts with a term which vanishes at the 
(12) cluster-separated limit. The (12)-limit equation is then built with the ap- 
proached — 2i7r/?i/32ifi2 only. The cluster separability can however be restored 
by performing a supplementary iteration, like our repl acement of V by Pq—S in 
(P7), leading the the cluster-separable correction (|l03| ). Another possibility con- 
sists in our previous choice (|6l| ) of . With this choice, the cluster-separated 
limits are exact. Moreover, the corresponding two-body scattering amplitudes 
are also correct. Our choice (^) of is thus a priori the most interesting. A 
simpler choice (as its Born approximation) would however not be uncorrect in 
the calculation of the higher-order corrections to the three-fermion bound state 
spectrum. 



4 Conclusions 

The kernel-approximated reduction is an interesting alternative to the more 
usual propagator-approximated reduction, although the resulting 3D equation 
is valid only for bound states. In both cases, the 3D potentials can be ob- 
tained with simple modifications of the Feynman rules which give the transition 
operator in field theory. Further investigations on the kernel-approximated re- 
ductions, with comparison with the other method (one-body limit, introduction 
of an external potential, gauge invariance, numerical examples, etc..) would be 
useful. 
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Our 3D reduction of the three-fermion Bethe-Salpeter equation provides in 
principle a way of progressively improving a manageable approximated 3D equa- 
tion. Wc use it to give a manageable expression of the first-order energy shift. 
The contributions of the free propagator poles to the energy shift are easily 
calculated in terms of the Bethe-Salpeter kernels with two of the fermions on 
their mass shell, in the spirit of Gross' spectator model. The contributions of 
the singularities of the kernels themselves should be computed in more specific 
frameworks ( He^ atom, triton, three-quark systems...) and could be important, 
as suggested by our computation of the first-order energy shift in the framework 
of the one-photon exchange model for the two-fermion problem. 
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